One-loop amplitudes for four-point functions with two external massive quarks and 

two external massless partons up to 0{e^) 
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We present complete analytical 0{e^) results on the one- loop amplitudes relevant for the next- 
to-next-to-leading order (NNLO) quark-parton model description of the hadroproduction of heavy 
quarks as given by the so-called loop-by-loop contributions. All results of the perturbative cal- 
culation are given in the dimensional regularization scheme. These one-loop amplitudes can also 
be used as input in the determination of the corresponding NNLO cross sections for heavy flavor 
photoproduction, and in photon-photon reactions. 
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I. INTRODUCTION 

At the leading order (LO) Born term level, heavy quark 
hadroproduction has been studied some time ago Ij. 
The next-to-leading order (NLO) corrections to unpolar- 
ized heavy quark hadroproduction were first presented in 
0, 13 , and in 0, 13 for photoproduction. Corresponding 
results with initial particles being longitudinally polar- 
ized were calculated in |3 and tili|9|,ll0|. A calcula- 
tion of the NLO corrections to top-quark hadroproduc- 
tion with spin correlations of the final top quarks was 
performed in Analytical results for the so called 

"virtual plus soft" terms were presented in 0, 0, |3 for 
the photoproduction and unpolarized hadroproduction 
of heavy quarks. Complete analytic results for the po- 
larized and unpolarized photoproduction, including real 
bremsstrahlung, can be found in [Tol . 

It is well known that the NLO QCD predictions for the 
heavy quark production cross sections suffer from theo- 
retical errors because of the large uncertainty in choos- 
ing the renormalization and factorization scales. In spite 
of considerable progress due to recent work in bringing 
closer theory and experiment (see e.g. 0, the 
need for next-to-next-to-leading order (NNLO) results 
for heavy quark production in QCD is by now clearly 
understood. The NNLO corrections are expected to sig- 
nificantly reduce the renormalization and factorization 
scale dependence inherent to the NLO parton model pre- 
dictions. 

During the last several years much progress has been 
achieved in developing and applying various techniques 
for an all order resummation of heavy quark production 
cross sections in different reactions. This concerns the 
resummation of the divergent terms in some specific re- 
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gions of phase space (so called large logarithms) to NLO 
(NLL logs) and NNLO (NNLL logs) leading logarithmic 
accuracy. We may mention the work on the threshold 
and recoil resummations of NLL logs in hadronic 
collisions. Much activity was also devoted to the resum- 
mation of NNLL threshold logs for heavy quark produc- 
tion in e^e~ (see e.g. the informative review '^15] and 
references therein) and 77 reactions. However, this 
cannot replace the need of having the exact NNLO re- 
sults for obvious reasons. In fact, these resummed results 
could be better understood when the exact NNLO results 
are available. 

The full calculation of the NNLO corrections to heavy 
hadron production at hadron colliders will be a very diffi- 
cult task to complete. It involves the calculation of many 
Feynman diagrams of many different topologies. It is 
clear that an undertaking of this dimension will have to 
involve the efforts of many theorists. As one example, 
take the recent two-loop calculation of the heavy quark 
vertex form factor ^3 which can be taken as one of the 
building blocks of the NNLO calculation. Another build- 
ing block are the so-called NNLO loop-by-loop contribu- 
tions which we have begun to calculate. The necessary 
0{e^) one-loop scalar master integrals that enter the cal- 
culation have been determined by us in . The present 
paper is devoted to the determination of the correspond- 
ing ©(e^) gluon- and quark-induced one- loop amplitudes 
including the full spin and color content of the problem. 
In a sequel to this paper we shall present results on the 
square of the one-loop amplitudes thereby completing the 
calculation of the loop-by-loop part needed for the de- 
scription of NNLO heavy hadron production. 

In Fig. n we show one generic diagram each for the four 
classes of contributions that need to be calculated for the 
NNLO corrections to the gluon-initiated hadroproduc- 
tion of heavy flavors. They involve the two- loop contri- 
bution 1^), the loop- by-loop contribution 1^), the one- 
loop gluon emission contribution (^) and, finally, the 
two gluon emission contribution (^i). The correspond- 
ing graphs for the quark-initiated processes are not dis- 
played. 



a) 



b) 




c) d) 
FIG. 1: Exemplary gluon fusion diagrams for the NNLO calculation of heavy hadron production. 



In this paper we concentrate on the loop-by-loop con- 
tributions exemplified by Fig.^. Specifically, working in 
the framework of the dimensional regularization scheme 
[igi. we shall present 0{e^) results on the one-loop am- 
plitudes. The expansion of the one-loop amplitudes up 
to £^ is needed because the one-loop integrals exhibit ul- 
traviolet (UV) and infrared (IR)/collinear (or mass(M)) 
singularities up to 0(e^^). When squaring the one-loop 
amplitudes to obtain the singular and finite parts of the 
loop-by-loop contributions one must thus know the one- 
loop amplitudes up to e^. 

In dimensional regularization there are three different 
sources that can contribute positive e-powers to the Lau- 
rent series of the one-loop amplitudes. First, one has 
the Laurent series expansion of the scalar one-loop in- 
tegrals which have been calculated up to 0{e^) in [l8l| . 
Second, the evaluation of the spin algebra of the loop am- 
plitudes brings in the n-dimensional metric contraction 
9fj.i^9^'^ — n — 4 — 2e. Third and last, the Passarino- 
Veltman decomposition of tensor integrals will again 
bring in the metric contraction g^vQ^'^ = n = 4 — 2e. 
The latter two points will be treated in this paper. It 
is clear that through the interplay of the three differ- 
ent sources of positive e-powers the Laurent series of the 
one-loop amplitude itself will, order by order, contain 
different orders of the Laurent series coefficient of the 
scalar integrals. 

We have confirmed the results on the Laurent expan- 
sion of the one-loop amplitude up to 0{e^) presented in 
p^. These results will not be listed again in this paper. 
In this paper we present analytical results for the coeffi- 
cients of the £- and e^-terms of the e-expansion includ- 
ing also their imaginary parts. When presenting our re- 
sults, we shall make use of our notation for the coefhcient 
functions of the relevant scalar integrals calculated up to 
0{e'^) in [13. For the calculation of the one-loop dia- 
grams with two external massive quarks and two external 
massless partons one needs one scalar one-point function 
A, five scalar two-point functions -Bi, six scalar three- 
point functions Ci, and three scalar four-point functions 
Di. For example, for the scalar four-point functions Di 
we defined successive coefficient functions Df according 



to the expansion 

+£2Z?f' +0(£3)}, (1.1) 

where Ce{m'^) is defined by 

Similar expansions hold for the scalar one-point function 
A, the scalar two-point functions Bi and the scalar three- 
point functions Ci. For the convenience of the reader we 
have included a table from jl^l where all the necessary 
one-loop master scalar integrals are listed. We note that 
for the one-loop scalar integrals the UV and IR/M sin- 
gularities never overlap, i.e. do not multiply each other. 
Singularities of order appear only when both IR and 
M poles are present simultaneously. This last case is real- 
ized when the massless gluon is attached to either mass- 
less fermion or a gluon line in the Feynman diagrams. 
Consequently, graphs ^1), ©l), (Efl) and (|3f2) shown 
in the next section have only poles, while graphs 
(jSli2), |(3Ji3), ijStS) and ^2) have e"^ poles. The details 
of the pole structure of the various Feynman diagrams 
can be found in 203 . 

As remarked on before we have endeavoured to calcu- 
late the loop-by-loop contributions in three steps starting 
with the scalar one-loop integrals, then calculating the 
one-loop amplitudes and finally squaring the one-loop 
amplitudes. If one's interest is only in the unpolarized 
rate one can directly move from step 1 to step 3 with- 
out the interim step of having to evaluate the one-loop 
amplitudes. However, in the latter case one loses the 
information on the spin content of the one-loop contri- 
butions which cannot be reconstructed from the rate ex- 
pressions. On the other hand, having expressions for the 
one-loop amplitudes allows one to easily derive the one- 
loop contributions to partonic cross section including any 
polarization of the incoming or outgoing particles. Our 
results on the one-loop amplitudes are given separately 
for every Feynman diagram in order to facilitate the use 
of the results for other relevant processes that differ by 
color factors. 
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TABLE L List of one-, two-, three- and four-point massive one-loop functions calculated in our previous paper [T^ up to 0{e^). 





Nomenclature of f3] 


Our nomenclature Novelty Comments 


1-point 


A{m) 


A 


Re 


2-point 


B{pi -p2,0,m) 


Bi 


Re 




B{pz +p4,m,m) 


B2 


- Re, Im 




B{p4, 0, m) 


Bs 


Re 




B{p2, m, m) 


Bi 


Re 




B(P3+P4,0,0) 


Bs 


- Re, Im 


3-point 


C(p4,P3,0,m,0) 


Ci 


new Re, Im 




C(p4, ~p2,0,m,m) 


C2 


new Re 




C(— p2,P4,0, 0,m) 


Cs 


Re 




C(-P2, -Pi, 0,0,0) 


C4 


- Re, Im 




C(-p2, ~pi,m,m,m) 


Cs 


- Re, Im 




C(p3,P4,m,0, m) 


Ce 


- Re, Im 


4-point 


D{p4, -p2, —pi, 0, m, m, m) 


Di 


new Re, Im 




D(-p2,P4,P3,0,0,m,0) 


D2 


new Re, Im 




D{-p2,p4., -pi,0,0,m,m) 


Di 


new Re 



The hadroproduction of heavy flavors proceeds 
through the following two partonic channels: 

g + 9^Q + Q, (1-3) 

where g denotes a gluon and Q{Q) denotes a heavy quark 
(antiquark), and 

q + q^Q + Q, (1.4) 

where q[q) is a light massless quark (antiquark). 

Note that the Abelian part of the NLO result for (|1.3() 
provides the NLO corrections to heavy flavor production 
by two on-shell photons 

1 + + (1.5) 

with the appropriate color factor substitutions. The re- 
sults for (|1.3() can also be used to determine the corre- 
sponding amplitudes for heavy flavor photoproduction 

1 + + (1.6) 

We mention that the partonic processes (|1.3|l and (|1.4|l 
are needed for the calculation of the contributions of 
single- and double-resolved photons in the photonic pro- 
cesses (|1.5|l and (|1.6|l . 

NLO cross sections for the process (|1.5f) have been de- 
termined in m mm for unpolarized and in [H 111] 
for polarized initial photons. Note that the authors of 
|24| used a nondimensional regularization scheme to reg- 
ularize the poles of divergent integrals. In the papers 
pH 01 analytic results were presented for "virtual plus 
soft" contributions alone. We also note that complete an- 
alytical results including hard gluon contributions can be 
found only in . The two-photon reaction (|1.5|) will be 
investigated at future linear colliders. NLO corrections 
for the heavy quark production cross section (|1.5|) with 



incident on-shell photons in definite helicity states are 
of interest in themselves as they represent an irreducible 
background to the intermediate Higgs boson searches for 
Higgs masses in the range of 90 to 160 GeV (see e.g. 
|23ll24| and references therein). 

The paper is organized as follows. Section ITU contains 
an outline of our general approach as well as one-loop 
amplitudes for the gluon fusion subproccss for the self- 
energy and vertex contributions including their renormal- 
ization. In Section III we discuss the one-loop contribu- 
tions to the four box diagrams in the same gluon-gluon 
subprocess and give a detailed description of our global 
checks on gauge invariance for our results. Section IV 
presents analytic results on the quark-antiquark subpro- 
cess (11.411 . Our main results are summarized in Section V. 
Finally, in two appendices we present results for the var- 
ious coefficient functions that appear in the main text. 



II. CONTRIBUTIONS OF THE TWO- AND 
THREE-POINT FUNCTIONS TO GLUON 
FUSION 

The Born and the one-loop contributions to the par- 
tonic gluon fusion reaction (7(^1) -I- 3(^2) — > Q{p?,) + Q{pa) 
are shown in Figs. In this section we discuss our 

evaluation of the self-energy and vertex graphs that con- 
tribute to the above subprocess. With the 4-momenta 
Pi{i — 1, 4) as shown in Fig. Inland with m the heavy 
quark mass we define: 

s=(pl^-p2)^ t = T -n? = {pi-p'if -m?, 

u = U — m? = (j)2— PsY — m?- (2.1) 

In order to isolate ultraviolet (UV) and in- 
frared/coUinear (IR/M) divergences we have carried out 



4 



Pi \l 

b 

a 

P2 V 



P3 



Q 



P4 



Pi 



P2 



Pi 



P2 




FIG. 2: The f-, u- and s-shannel leading order (Born) graphs contributing to the gluon (curly lines) fusion amplitude. The 
thick solid lines correspond to the heavy quarks. 



all our calculations in the dimensional regularization 
scheme (DREG) 19] with the dimension of space-time 
being formally n — A ~ 2e. 

First of all we note that in general the amplitudes for 
all the Feynman diagrams in the gluon fusion subprocess 
can be written in the form 

M = eMMP2MP3)M^''v{p4), (2.2) 

For purposes of brevity, we will present our results in 
terms of truncated amplitudes M^'^ where the polariza- 
tion vectors and Dirac spinors are omitted. For reasons 
of brevity we shall also refer to the truncated amplitudes 
as amplitudes. Of course, the presence of the polariza- 
tion vectors and Dirac spinors is implicitly understood 
throughout this paper in that the mass shell conditions 
Pi^fiiPi) = and ^3u(p3) = muij)^) etc. are being used 
to simphfy M'^'^ [2£]. Furthermore, M^"^ contains the 
common factor C^iw?) defined in Eq. \i.2l which arises 
from the scalar one-loop integrations described in [l8l| . 
Throughout the paper we will omit from all our one-loop 
M^'' amplitudes the common factor 

C = g^C,{m'), (2.3) 

where g is the renormalized coupling constant. 

There are three sets of contributing graphs: The t- 
channel, u-channel and the s-channel graphs as exempli- 
fied in Fig. 121 for the LO Born term contributions. Since 
the M-channel amplitudes can be obtained from the 
f-channel amplitudes Mt by the relation 

Mt'^Mu^ia^h, pi^p2, fJ,<^iy}, (2.4) 

we shall not list results of the w-channel contributions. 
In 1)2. 4|l a, b are the color indices of the two gluons. We 
make it clear from the outset that additional u-channel 
graphs are obtained from the relevant ^-channel graphs 
by the interchange of the two external bosonic lines (not 
only momenta). In exception are the two vertex inser- 
tion diagrams (3c3) and (3c4) which will be discussed 
later on. All three interchanges (color, Lorentz indices 
and bosonic momenta) have to be done simultaneously. 
Note that the second interchange in H2.4|l implies also the 
interchange t <^ u. In general, when speaking about the 
t-u symmetry of a given subset of amplitudes, we will 
imply invariance of those amplitudes under the transfor- 
mations (E3I). 



We start by writing down amplitudes for the leading 
order Born terms. For the ^-channel gluon fusion subpro- 
cess (first graph in Fig. 12) we have: 

Br = ^iT'T^riJ^s -h + m)77i, 

where and T" are generators (T" = /2, a = 1, 8 
and the A° are the usual Gell-Mann matrices) that de- 
fine the fundamental representation of the Lie algebra of 
the color SU(3) group. Analogously, for the u- and s- 
channels depicted in the second and third graph of Fig. [21 
we have, respectively, 

where the tensor C[^'^" is obtained from the Feynman 
rules for the three-gluon coupling and is given by 

Ca""^ = 9i^u{jPi ~P2)a+9va{pi + 2p2)m -3m'^(2pi +P2)u- 

(2.5) 

We have omitted a common factor in the Born am- 
plitudes. Acting with Dirac spinors u{jp^) and v{p4) on 
the above truncated Born amplitudes from the left and 
the right, respectively, and using the effective relations 
Pi = P2 ~ 0, as remarked on before, we arrive at the 
following expressions for the leading order amplitudes: 

Br = iT'T^r^iY - "^P^sll/t; 

B^ = 2i{T''T'' -T''T''){g'"'^i+pt^j'' -p'(J^')/s. 

Next we proceed with the description of the two-point 
insertions to the amplitudes of the subprocess 1)1. 3|) . But 
before we turn to the two-point functions one should 
mention that our choice of renormalization scheme will be 
a fixed flavor scheme throughout this paper. This implies 
that we have a total number of flavors rif = nif + 1, where 
nif is the number of light (i.e. massless) flavors and the 
"1" stands for the produced heavy flavor. Thus there will 
only be nif light flavors involved/active in the (3 function 
for the running a QCD coupling as, and in the splitting 
functions that determine the evolution of the structure 
functions. When having massless particles in the loops 
we are using the standard MS scheme, while the con- 
tribution of a heavy quark loop in the gluon self-energy 
with on-shell external legs is subtracted out entirely. 
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FIG. 3: The t-channel one-loop graphs contributing to the gluon fusion ampUtude. Loops with dotted hnes represent gluon, 
ghost and light and heavy quarks. 



Consider first the two f-channel self-energy insertion 
graphs (3d2) and (3d3) in Fig. 3 with external legs on- 
shell. These graphs are very important as they deter- 
mine the renormalization parameters in the quark sec- 
tor. Throughout this paper wc use the so called on-shcU 
prescription for the renormalization of heavy quarks, the 
essential ingredients of which we describe in the follow- 
ing. When dealing with massive quarks one has to choose 
a parameter to which one renormalizes the heavy quark 



mass. It is natural to choose a quark pole mass for 
such a parameter - the only "stable" mass parameter 
in QCD. The condition on the renormalized heavy quark 
self-energy Sr(j?') is 

Sr(^)|j/=m = 0, (2.6) 

which removes the singular internal propagator in these 
self-energy insertion diagrams. This can be seen from the 
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FIG. 4: The s-channel one-loop graphs contributing to the gluon fusion amphtude. Loops with dotted hnes as in gl,h,jl and 
j2 represent gluon, ghost and light and heavy quarks. The four-gluon coupling contribution appears in g2. 



explicit result for the renormalized heavy quark exter- 
nal self-energy Er(i') e.g. in dimensional regularization 
scheme: 



1 + 



£(1 - 2e) 

2 ^ 
m — p 



2p2 



-(!-£) 



? 2 

p — m+ p X 



9 9 
m (2 — e) 



(2.7) 



The above condition H2.6() determines the mass renormal- 
ization constant Zm- For the wave function renormaliza- 
tion we have used the usual condition (see e.g. Ref. 



d 



(2.8) 



which fully determines the wave function renormalization 
constant Z2. Since the condition 1)2. 8|l is not mandatory 
in general, there is a freedom in determining the constant 
Z2. Note that the condition 1)2. 8|l sets all external heavy 
quark self-energy insertion diagrams to zero, thus making 
the heavy quark case similar to the massless one in this 
regard. Below we list our expressions for the mass and 
wave function renormalization constants. 
In the DREG scheme we arrive at the result 



1 - g^CFCim^) 



2e 



e(l - 2e) 



(2.9) 



which can be expanded in e to give 



Z2 — 'Z'rai 



(2.10) 



where (7^^=4/3 and we do not make a distinction which 
poles are of ultraviolet or IR/M origin as we did in |20|. 
After the mass renormalization procedure is applied we 
obtain the final results for the two self-energy insertion 
graphs in the DREG scheme 



'(3d2) 



'(3d3) 



3-2e 
e(l - 2e) 



(2.11) 



-CF^r ( ^ + 4 + 8e + IQe^ + 0{e-^ 



From here on we will present only results for the e and 

order contributions to the amplitudes. 

After addition of the mass renormalization countert- 
erm the contribution of the quark self-energy insertion 
graph (3dl) with external legs off-shell reads: 



fc=i 



+Bf-^h/T - kl&m^lt) 
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k=l 



-B 



(^-^^ /T - k8/tj 



(2.12) 



The coefficients b['^'^ and B\^"'^' come from the Lambent 
series expansion of the scalar two-point function Bi (see 
Table 1) quite similar to the corresponding Laurent se- 
ries expansion of the four-point functions Di shown in 
Eq. lfLT|l . 

The remaining quark self-energy insertion graphs (4il) 
and (4i2) with external on-shell legs are derived in anal- 
ogy to the ones considered above: 



^"(4il) 



^"(4i2) 



'CfB^" (8e + 166^ 



(2.13) 



Concerning the gluon self-energy insertion graphs (3el) 
and (3e2) with external legs on-shell, the only nonvanish- 
ing contributions are those from heavy quark loops. They 
are given by 



(3el) 



1 2 
e 3' 



(2.14) 



However, these contributions are explicitly subtracted 
(together with the common factor C^{m?), see Eqs. (|1.2f) 
and H2.3|l ) in the on-shell renormalization prescription. 
Therefore, due to the UV counterterm that subtracts 
this loop with heavy quarks, there are no finite contribu- 
tions to the amplitudes from these self-energy diagrams. 
However, at the same time this counterterm introduces 
the pole terms from the light quark loop sector that are 
needed to cancel soft and coUinear poles from the other 
parts of the amplitude, e.g. from the real bremsstrahlung 
part. This indicates that in practice it is very hard to 
completely disentangle UV and IR/M poles in heavy fia- 
vor production and in most cases one obtains a mixture 
of both instead. 

For the reasons specified above we present the gauge 
field renormalization constant Z3, used for the gluon self- 
energy subtraction: 



1 



e 



= l + -^^{Po- 2Nc)C,{fi') - -a(m^)| , (2.15) 

where the QCD beta-function Pq — {llNc — 2nif)/3 con- 
tains only light quarks. Nq = 3 is the number of colors. 
Accordingly, for the coupling contant renormalization we 
obtain 



1 



e 



f3o 



-C,(m2 



(2.16) 



As was the case for the diagrams (3el) and (3e2), di- 
agrams (4jl) and (4j2) also vanish altogether due to the 
explicit decoupling of the heavy quarks in our subtrac- 
tion prescription. However, instead of renormalizing sep- 
arately each Feynman diagram, one can chose to employ 
the renormalization group invariance of the cross section 



and do only a mass and coupling constant renormaliza- 
tion. In this case, knowing the results for the gluon self- 
energy diagrams turns out to be useful in checking the 
complete cancellation of UV poles by just rescaling the 
coupling constant in the LO terms f/barc — * Zgg. One has 



,1 2 
e 3' 



(2.17) 



Finally we arrive at the gluon self-energy insertion 
graph (4h) , which contains the off-shell gluon self-energy 
loop that is used for the derivation of the renormalization 
constant Z3. We have evaluated the internal loop in the 
Feynman gauge. In our result we show separately the 
gauge invariant pieces for gluon plus ghost, light quarks 
and one heavy quark flow inside the loop: 



Bi" 



1 



]\T n — 14 



-8e 



2e) 



2(l-£) 



3-2e 



1 2 



(2.18) 



with n = 4— 2e in the DREG scheme. Be, is the two-point 
integral whose explicit form is given in Pt- We expand 
the first line of H2.18|l in powers of e and find 



Bi^" 



c 



1 5 
e 3 
1132 31 

y 

10 



31 



188 
"27 



■C(2) 



10 



81 
1 2 
e 3 
328 



C(2)-^C(3) 



56 
27 



■C(2) 



10 4 

81 -T^(^)-3^(^) 



— S 
m 2 



1 2 

£ 3 



with 
/ = 1 



B. 



(0)3-^ 



(2.19) 



(2.20) 



3 

15(0)^2^^(1) 3^ 
932^2 2 ^ 

.±_2s(");32_ 15(1)^2^5(2) 3-^'^ 
27 9 2 3 2 2 2 



In (|2.20|l we have made use of the definition 



13 = y^l - Am?/s. 



(2.21) 



Concluding our discussion on the 2-point insertions 
we remark that the amplitudes for the relevant u- 
channel 2-point insertion diagrams can be obtained from 
Eqs. l|TTT|l . and lfni|l by the transformation ((231 • 

Next we discuss the t- and w-channel vertex insertions. 
In this paper we write down only the e- and £^-terms of 
the Laurent expansion. The terms proportional to £^^, 
and £° can be found in [23|- We begin with the purely 
nonabelian graph (3b) with the four-gluon vertex. The 
amplitude takes the following form 
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^^-'(3b) 



= iNc{T''T' e''{(2P37^ + pIi"" - vh" - '^V%l"){Bf^ + 2Cf ^m' - 4fc) - mg^'"' {2Bf^ + 25, 
fc=i 

+4Cf ^m^ + Cf "^\s - 12fc) + 3m7^7'^(2Bf ^ + cf - 8fc)/2}/(s/?2) + (a <-> &, ^ ^ ;/)) 

2 

+ ^<5'^''^e'={(p!^7^ -p:^7^ +P^7'' + 2Cf 'm^ - 4fc)/2 



(fc-i) 

5 



fc=l 



I 



(2.22) 



It is easily seen from Eq. H2.22|l that the amphtude for 
the graph (3b) is exphcitly t-u symmetric, as it fohows 
from the geometric topology of this graph. It is thus 
important to state that there is no u-channel equivalent 
of graph (3b). 

Next we turn to graphs (3cl) and (3c2). As mentioned 
before, these diagrams occur also in other processes such 
as photoproduction and 77 production of heavy flavors 
when one or two of the gluons are replaced by pho- 
tons. For this reason we also present the correspond- 
ing t-channel color factors for these graphs. Then it is 



straightforward to separate our Dirac structure from the 
color coefficients and one can easily deduce the corre- 
sponding results for the other processes involving pho- 
tons. In order to facilitate this transscription we list the 
color factor for both diagrams (3cl) and (3c2) which turn 
out to be the same: 



T, 



(3cl) _ y(3c2) 



-A 



{Cf - 11£)t''T'' = -IT'T". (2.23) 



1 



col 



6 



The complete amplitudes are: 



M(7^^) = £''{B^i\'^Tn'/t + 1) + 2B^^-^\.t/t + 2&^^m' + 467^"'^™^ _ k%{Am'/t + l)}/6 

fc=i 

2 

+ iT^T'^ip'^Y ^^{Bi^ i^'t/t + t/T) + {2zt/t - t/T) + 2(cf ^ + 2Cf V - kSzt/t} 

k=l 

2 

+ mp^Al" ^^{Bi^lT - Bf'^\2/t + l/T) - 2Cf + M/t} 
fe=i 



m^^^^- J2 £'={Bf' + Bp'' + "''i - 6fc})/(30, 



fc=i 



where we have introduced the abbreviation zt = 2m? ~\-t. For the graph (3c2) we obtain: 



(2.24) 



^(3c2) = Sr^£''{Sr'(6mV< + l) + 2Bf"''ztA + 2Cf'm2-|-4C^ 
fe=i 

2 

+ tT''T''{plj'' J2 e'^lsf H-2mVi - 3 + t/T) - B[''-^k/T - 20. 



k8{Am'^/t+l)}/6 



(fc) 2 



k8T/t} 



fc=i 



+ mp2(2p^ - 7'^A) ^e'={i?fVT^ - B['"'\2/t + l/T) - 2Cl 



(/c-l) 



M/t} 



^7^7" '"'^^^i^ + ^^''"'^ + ~ 6fc})/(3i). 



(2.25) 



fc=i 
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Next we write down the results for graphs (3c3) and 
(3c4). The color factors for both diagrams are the same: 



(3c3) 



We have 



T, 



(3c4) 



col 



_r-phrjia 

2 



(2.26) 



fe=i 



+M(3mVi + 1)} 

2 

+3ir^r°(p^7'' ^ £'={B|''^m2(l/T - 2/i) 

k=l 

+B\^~^h/T - C^^h + kAzt/t} 

2 

+3m7^7" ^ £'={Bf ^ /2 - 2fc} 



2 



+mp^,i>iYY.^''{B^^\2/t-l/T) 
+s(''"^V7'-fc8/i})A. (2.27) 



And 



^(3c4) = 3Br^e'={-3Bf)mVi-Q'''t 



/c=l 



-fc4(3mVi + 1)} 

2 



k=l 



(fc-i). 



'(t/T - 2) + + fc4(2mVt - 1)} 

2 

+3m7''7'^ ^ £'={bJ'=V2 - 2fc} 
fc=i 

2 



fc=i 



(2.28) 



The results for the amplitudes of the relevant u- channel 
vertex insertion diagrams are obtained from Eqs. H2.24|l . 
(103), and by the transformation ^T^. 

However, there is a subtle point involved here: we stress 
that for the graphs (3c3) and (3c4) the Mt <-> M„ trans- 
formation (|2.4|) transforms the t-channel result of the 
graph (3c3) to the u-channel result for the graph (3c4), 
while the i-channel result of (3c4) goes to the u-channel 
result for (3c3). This is important to keep in mind when 
dealing with reactions which involve asymmetric set of 
graphs as e.g. in the photoproduction of heavy flavors. 
The reason for this is that when doing transformation 
H2.4(l the three-gluon vertex attached to one of the initial 
bosonic lines does not stay attached to the same bosonic 
line. However, we note that transformation's <-> p4 does 
uniquely relate all the t- and u-channel diagrams for the 
subprocess under consideration. 



Next we turn to the remaining s-channel graphs shown 
in Fig.^ For all the gluon propagators we work in Feyn- 
man gauge. This set of graphs is purely nonabelian for 
the QCD type one-loop corrections. In the case that one 
wants to replace the gluonic vertex correction in graph 
(4fl) by a photonic vertex correction one needs the ex- 
plicit form of the color factor for graph (4fl): 



(4fl) 



Nr. 



1 



col 



The amplitude including the color factor is 



(2.29) 



^'^(4fl) 



E £'"{3^2''^ + 2Bf + ^s(l + P^) 
fc=i 

-16/c}/6 -I- 2i{T''T'' - T''T'')m[-g'"'{s + 2t) 
2 

k=l 

/{&s^(3^). (2.30) 



Graph (4f2) contributes as: 



^(4f2) = A^cSrE^'i^s + 
fc=l 

-fcl6(5m^ - s)}/{2sl3^) 
+2iNc{T''T^ - T''T'')m[~gt"'{s + 2t) - 
2 

fc=i 

— 2£Sg s + K>Ll m s — CI s — 
k4{12m^ ~ s)}/2s^(3'^. (2.31) 

We end our consideration of the vertex insertions for 
gluonic fusion with the sum of the two graphs (4gl) and 
(4g2) which we refer to as the triangle graph contribution 
(tri)=(4gl)-|-(4g2). For the case when one has gluons and 
ghosts inside the triangle loop we obtain: 



M^^-i)(9) = -37Vc(i?rE^'i207BrVl2B 



j(fc-i) 



fc=i 



+54Cf + 8k+(k- 1)85^°^} 
2 

_^g-(yay6 _ T^T")^! E e''{g^"'[9B, 



(k) 



k=l 



-l2Bf^^^ + 9Cf\s -8k-{k- l)8B^°Vs 

+8p^K[3Bf + 2k + {k- l)2B^°Vs'} 
)/324, (2.32) 

where B^^^ = B^'' — 4/3. When one has light and heavy 
quarks inside the loop one has 

2 



,(fc-l) 



2k 



k=l 
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-{k - 1)2B^°^} - ii{T''T^ - T'T")^! X where nif is the number of Ught flavors in the triangle 
2 loop. For the heavy flavor case one has 

fc=i 

+ (fc-l)(5B^"^+3)})/81 (2.33) 
I 

2 

^rtri)(0) = 6(Sr ^^{^(^B^2^ + 2Bf "'^ + (fc - l)4S^"V3)mVs + ^ - 3sf - 2fc - 2(fc - 1)5^"^} 
fe=i 

2 

_^^yay6 _ T^'T'^) A [g^7s - 2p^K/s^] e''{12(3Bf ^ + 2sf + (/c - l)4S^"V3)mVs 

fc=i 

+18(cf ^ + Cf + (fc - l)cf V + 35^"'^ + 5fc + (A: - 1)(5S^''^ + 3)})/81, (2.34) 



where i?2 =^24/3. The complete amplitude for 
the triangle (tri)=(4gl)+(4g2) is the sum of the above 
three expressions (IT^ and ifT^ : 

^rtH) = ^rtH)(5) + A^rtH)!?) + Af^'^^ijiQ)- (2.35) 

In Ref. j2^ one can find general results for the gluon 
triangle in any gauge and dimension. We have compared 
the first two terms in H2.35|l with the corresponding ex- 
pressions in Ref. |2^ and found complete agreement. 

III. RESULTS FOR THE BOX DIAGRAMS IN 
GLUON FUSION 

In this section we describe the technically most in- 
volved derivation of the 4-point massive box diagrams. 
The four box graphs (3al)-(3a4) contributing to the sub- 
process g + g ^ Q + Q are depicted in Fig. 13 We have 
used Passarino-Veltman techniques to reduce tensor 
integrals to scalar ones where the scalar master integrals 
are taken from our previous publication . 



For each of the gluon fusion box diagrams we expand 
the truncated amplitude Af^" in terms of a set of 20 
Lorentz-Dirac covariants multiplied by the same number 
of invariant functions. In the reduction of the Lorentz- 
Dirac structure to this basic set of 20 covariants we have 
been making use of the mass shell conditions described 
in Sec. II. The 20 Lorentz-Dirac covariants are subdi- 
vided into eight subsets according to their Dirac struc- 
ture. The 20 invariant functions multiplying the covari- 
ants are sorted according to the contributions of a basic 
set of functions /f^' (called basis functions) related to 
the scalar master integrals of |18| . The index i runs over 
the members of the set of basis functions occuring in a 
particular graph. The index k denotes the power of e 
which the basis function multiplies. The basis functions 

(k) (i) 

are multiplied by coefficient functions h^J^ where the 
index pair (n, j) identifies the covariant which the coef- 
ficient function multiplies. Note that the basis functions 
/j*-*^-* have been defined such that the coefficient functions 

&)„ do not depend on the index k. We thus cast the box 
amplitude into the following universal form: 



M'^- = ^r™i Ee^ACE/'^^'?^ (3-1) 

+ A[5-E/^'^^^i^ +^^T.f't'^ +p^pij:f^^ +p^P3 Ei^'^^ii^ +p^pij:f''^^^ 

+ 7"b3 E ft'^ +p'^T. f^^'b^] + 7>!^ E f^^'b^ +p'^i: ft'b^] + -^rr E ft'^ 

+ [pS E f^""^ + P4 E f^"'b^^ + K E f^"'b^ +P'^i: 

+ m[5-Ei^"C^ +p^p^Y.f'''^ +p^i^f^^ +p^^Y.f't'^ +p^p^f?^^w 

+ {Mt ^ Mu}- 



The symbol {Mt ^ -Mu} at the end of Ea. H3.1|) needs to be explained. It has the same meaning as the syni- 
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bol Mt Mu defined in Ea. H2.4|l except that dia- 
grams (3a3) and (3a4) are exempted from tfie sum. The 
crossed boxes (j2t3) and (|3i4) go into each other under 
the M.t ^ J^u operation. More exactly, for each of these 
diagrams, when one symmetrically interchanges the two 
bosonic lines (together with the appended three-gluon 
vertex) one arrives at the original box graph topology 
since these boxes represent diagrams of the so called non- 
planar topology. This becomes even more clear when one 
interchanges ^ P4: In this case each of the two crossed 
box graphs is reflected into itself. 

Taking parity into account one has altogether 2 • 2 • 
2-2/2 = 8 independent amplitudes and thus eight in- 
dependent covariants for the process g + g^Q + Qin 
n = 4-dimensions. We have made no attempt to reduce 
the 20 (plus 3 from the w-channel contributions) covari- 
ants in H3.1|l to a basic set of independent gauge invariant 
covariants. In fact, gauge invariance will be checked later 
on in terms of the expansion (|3.1|) . At any rate, the num- 
ber of independent gauge invariant covariants will very 
likely change going from n = 4 to a general n ^ A. 

Depending on the type of the box graph one has a 
different number of terms in the (i) summation in H3.1|l . 

(k) 

These numbers as well as the set of basis functions 
related to the scalar master integrals are specified below. 
The coefhcient functions b^j^ are given in Appendix A of 
this paper. 

In the expansion (|3.1() it is convenient to choose one 
covariant as the t-channel Born term amplitude struc- 
ture (and correspondingly a u-channel Born term 
amplitude structure). We define it as 



^3- h + -m^Y, 



(3.2) 



which, when taken between the spin wave functions im- 
plying the effective relations — 0, — 0; can be 
written as 



(3.3) 



For each of the box diagrams (P^l) and (Ol2) we found 
the following empirical relations between the and b'^^^ 
coefficient functions: 



d5) 



u(6) 



d5) 



(3.4) 



Because of the relations (|3.4f) we will not write down the 
results for the b[^ coefficients in the Appendix A. 

Next we present the color factors and basis functions 
for the abelian type box diagram |(21il). For this graph 
the sums over i in H3.1|l run from 1 to 17 for each of the 
20 terms. One has: 



cab 



Nr 



(3.5) 



(fe) 



a 



(fe-i) 



.(fc) _ ^(k) 

8 — '-'5 ' 



Ak) 
J9 



C, 



(fc-1) 



Ak) _ Mk) 

JlQ — ^6 ' 



Ak) _ p,(k-l) 
Jll — ' 



Ak) _ T^{k) 
J12 — ' 



Ak) 
J 13 



(3.6) 



(fe-2) 



= (fc- l)Df"2\ 

lU) 



The corresponding coefficient functions are listed 
in Appendix A. Many of the coefficient functions are in 
fact related to each other. One has 



dj) 



-tb 



U) 



-12« - -""lOn' .?V0. (3.7) 

And for any given values of n and j one has 



"11 n — ''"9n I 



6ll = .s/(2t)6ll, (3.8) 



foil = s(3yi2z,)b[t, 6i^„ = -stpyi2z,)b[t- 

Further relations are valid for particular sets of the 
parameters n, j: 



b\''l = s/{2t)bli, J =0,1,2,4; 



dj) 



b^;i^.s/{2t)b<^, 



n — 1,3, 5; 



(3.9) 



bi'^ = s/{2t)b'i:, j=4,5,6,7. 



Because of these relations among the coefficient func- 
tions we will write down only the independent coefficients 
b'fj^ in Appendix A. 

For the nonabelian box diagram ||3^2) the sums over i 
in l|3.1|l again run from 1 to 17 for each of the 20 terms 
in H3.1|l . For the color factor we obtain: 



(3.10) 



The relevant seventeen basis functions that describe the 
result of evaluating the box diagram ^^2) are given by 







Ak) _ 
J2 " 


- B^^'> 

- , 


Ak) 

J3 " 


Ak) 
Ji - 


R(fe) 


Ak) 
J5 — 


(jik~i) 


Ak) 
J& 


Ak) _ 
/? — 




Ak) _ 

Js - 


- ^3 ' 


Ak) _ 
J9 - 


Ak) _ 

Jw - 




Ak) _ 

Jll — 

Ak) 
Jl3 


^2 ' 
= fc. 


Ak) 
J12 



B 



(fc-i) 



CI 



(fe) 



ik~i) 



(3.11) 



f[l'> = {k-in 



ik-2) 



Ak) 
J16 



(fc-i)cf-'\ 



There are five relations between particular coefficients 
for the box diagram ||2^2), valid for any values of n and 

i- 

b^ = ./(2i)4^, b^£ = -(./2)4^, (3.12) 



and 



Ak) _ 

Jl ~ ^1 



p(fc) 



(fc) 



Ak) _ p,(fc) 

h - ^1 ^ 



Ak) 



c. 



ik-1) 



Ak) _ p,(fc-l) 



Ak) 
Je 



a 



(fc) 



= sz,/{2e)b[t, b^l = s/(2t)b^L 



b^^l = -is/2)b[l. 



(3.13) 
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In addition, one has two sets of relations that are vahd 
for the corresponding parts of the expression for the 
box ((3^2). The first set of relations is 



b'£^sz,/{21^)b^l 



"I5n ~ '^"7n ' 



(3.14) 
(3.15) 



■'l6n 



2b: 



9n ' 



2b 



(i) 



The above equalities are valid for j — 0,2,3,4,5,6 and 
for j = 1,7 and n = 1. Note that in the presence of 
the set (|3.13|l not all of the relations in (|3.14l) , H3.15|l are 
independent. Therefore, we can choose Eq. (|3.14l) and 
only one relation (e.g. the second one) out of the four 
relations in (|3.15|) as a set of independent relations. 

The second set of relations is represented by the two 
equalities that are identical to the ones of H3.12|l , but are 
valid only for j = 4, 5, 6, 7 or for j = 1 and n = 2, 3, 4, 5: 



bfl^s/{2t)b^i, 



b^l^-{s/2)bi^l 



(3.16) 



In the case of the crossed box (|3t4) one has twenty 
basis functions for each of the terms in 1)3.1(1 . The color 
factor for this graph takes the simple form 



col 



(3.17) 



The functions /f are defined as follows: 







Ak) 
J2 


= <, 


H 




Ak) 
Ji 


R(fe) 


Ak 
Jb 


= C^''"^' 


Ak) 

Je 


rik) 

— ^2 1 


Ak 

Jr 


Mk-l) 
— ^2u ' 


Ak) 
J8 


rik) 

— ^2u ' 


Ak 

Jd 


= c-c^-i) 


Ak) 
JlO 


rik) 

- ^3 ' 


Ak 

Jll 


Mk-l) 
— ^3u ' 


Ak) 
Jl2 


ri^) 


Ak) 
J13 


_ r,(k--i) 

— ^3 ' 


Ak) 
Jl4 


— ^3 ' 




Ak) 
il5 




(3 


[k- 




Ak) _ 
Jl7 - 


-{k-i)ct-^\ 


{k- 




Ak) _ 
J19 ~ 






fw — " 




2) 

7 



Ak) 
Jl6 

Ak) 
Jl8 



where the subscript "u" is an operational definition pre- 
scribing a {t ^ u) interchange in the argument of that 
function, i.e. sj^'' = B[''\t <^ u). 

There are numerous relations between the coeffi- 
cient functions for this diagram. These relations read: 
For any value of n and j 



"11 



dj) _ 

^13n — 



as well as 



di) 

^7n 



bi'^u/t, J + 5; 



(3.19) 



(3.20) 



42„ = 4W^, ^--4,5,6,7; 
&i2n = ^10 n"/*' n^\; 
b['L^-b[\\t, n^0,2; 

^(2) _ _^(2) 

Further one has a less general but still very useful relation 
for any n 

4n - -biiltu/{2D + tu), J - 0, 3, 6, (3.21) 

with D = nP's — tu. Equation (|3.21l) above is also valid 
for J = 1 , 7 and n — 1. 

For the coefficient functions that effectively only mul- 
tiply the £^-terms we have two sets of relations. One set 
is 



dj) _ 2hU) 

"16 n ~ ^"Sn ! 

Aj) _ ohii) Aj) _ ^dj) 

dl) 



(3.22) 



d3) do) _ Oh'-') 

11 n I '^20n — ^"13 n 



which are valid for the same values of j, n as specified in 
and after H3.21|l . The other set reads 

fell = b'fL^/t^ (3-23) 
b^^^^^-b^^l,tu/{2D + tu), 

fell = -feil"V(2i? + i"), 
4l = 4'„^"V(2^ + M- 

The relations (|3.23() are global for the crossed box (PJl4) , 
i.e. valid for any set of index values. Because the re- 
lations H3.22|) always occur together with the relations 
1)3.21(1 . only the first relation of 1(3.22(1 is important. The 
other four relations in 1(3.22(1 are redundant since they 
can be derived from 1(3.19(1 . the first relation in 13.20(1 . 

(1?:^ and in23- 

In addition to the relations listed above, various coeffi- 
cient functions of the crossed box are related by {t u) 
exchange. For instance, the coefficient functions mul- 
tiplying the Born term structure (or j = 0) are 
related by 



bf!!=b'f^2^,{t^u), i = 1,2, 5, 6, 9, 10; 
b'i!l=b[%,{t^u), z = 16,18; 



(3.24) 



The remaining (j — 0) coefficient functions turn into 
themselves under (i «-> u). 

Other coefficient functions are negatively related by 
(t <-> u)-exchange: 

C = -C2.3(^ ^ ^ = 1. 2, 5, 6, 9, 10; 



-b^^lit^u), i = 13,14,15,20; 
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FIG. 5: The lowest order Feynman diagram contributing to the subprocess qq QQ. The thick lines correspond to the heavy 
quarks. 



(3.25) 

^■"5^ - -^>2,2(i ^u), i = 1, 2, 5, 6, 9, 10; 
bf,l = -bfl{t^u), i = 13, 14, 15, 20; 

^■5 =-C\,2(i--"), ^-16,18. 

Furthermore, the whole term corresponding to = 4 in 
(|3.1|) is antisymmetric under {t <-> u). The foUowing 
pairs of coefficient functions are negatively related in the 
sense of (|3.25|l : the 6-2' are related to , and the 

are related to ^ where I can take any of the values 
I = i,i + l,« + 2 depending on the value of i. The number 
of independent coefficient functions is greatly reduced for 
this box because of all these relations. We took advantage 
of this fact when writing down the relevant coefficient 
functions in Appendix A. 

As explained after Eq. (|3.1I) the crossed box ((2^4) is 
obtained from (PJl3) with the help of the AAt <-> Mu 
operation. For this reason we write down explicit results 
only for one of the box ((3^3) in Appendix. A. 

A necessary check on the correctness of our one-loop 
results is gauge invariance. For example, for gluon 1 this 
implies that one must have 

Pi^,eu{p2)u{pz)M^''' [one - loop)i;(p4) = 0, (3.26) 

for each of the remaining independent amplitude struc- 
tures that multiply e.g. pii,, p^i, and 7^/. Similarly one 
must have 

P2.e^ipi)u{p3)M^''{one - loop)i>(p4) = (3.27) 

again, for each of the remaining independent amplitude 
structures that multiply e.g. p2/i, P3n and 7^. We have 
verified gauge invariance for the following gauge- invariant 
subsets of diagrams: (i) When the incoming gauge bosons 
are photons, i.e. including graphs ^^1), ijStl), l|2t2), 
(PJil), lI21i2), (|31i3) plus their u-channel counterparts 
with their corresponding color weights; (ii) For the photo- 
production of heavy flavors, i.e. including all the above 
diagrams plus graphs ^4), ^4), ^1) plus their u- 
channel counterparts, with corresponding color weights; 
(iii) For the hadroproduction of heavy flavors, which ul- 
timately includes all the graphs from Figs. Inland 0] plus 
their relevant u-channel counterparts. We emphasize 
that the above gauge invariance checks were made sepa- 
rately for both color structures Cp and Nc, and for every 



existing combination of color matrices T", and S"''', 
whenever they arise. When checking on gauge invariance 
all the relevant s-, t- and u-channel graphs have to be 
added. Gauge invariance must of course be checked for 
each power of e and for each of the coefficient functions 
of the Laurent series expansion of the scalar master in- 
tegrals separately, independent of their actual numerical 
values. 

Finally we note that the original computer output for 
the box diagrams was extremely long. The final results 
were cast into the above shorter form with_ the help of 
the REDUCE Computer Algebra System j23. 

IV. ANNIHILATION OF THE 
QUARK-ANTIQUARK PAIR 

The LO Born graphs contributing to this subprocess 
are shown in Fig. [5| In Fig. we show the graphs con- 
tributing at one-loop order. 

The leading order contribution proceeds only through 
the s-channel graph. One has: 

= iTt^T^MP2)ru{pi)u{p3Mpi)/s. (4.1) 

Here the color matrices T" belong to different fermion 
lines which are connected by the gluon having color in- 
dex a. We have again left out the factor in the Born 
term contribution (|4.1|l . In the Passarino-Veltman re- 
duction for tensor integrals we can make use of the same 
scalar integrals of ^3 as those appearing in the gluon 
fusion subprocess, with relevant shifts and interchanges 
of momenta when needed. 

Starting again with the 2-point insertions, we notice 
that the result for graph (6g) can be obtained from the 
one of (|2.19|) for graph (4h) in the gluon fusion subprocess 
by the simple replacement 

M(6g) = M(7,) (i?r - B,,). (4.2) 

The massless quark self-energy insertion graphs (6j) 
and (6k) with external legs on-shell vanish identically: 

Af(6j) = Af(6k) = 0. (4.3) 

The massive quark self-energy insertion graphs (6h) 
and (6i) with external legs on-shell are calculated in anal- 
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FIG. 6: The one-loop Feynman diagrams contributing to the subprocess qq — > QQ. The loop with dotted line represents gluon, 
ghost and light and heavy quarks. 



ogy to the ones considered in the previous section: 

3 — 2s 

M(6h) = M(6i) = -CpBgg 



s{l - 2e) ' 



The results for the vertex insertions are relatively 

short. Starting with graphs (6c) and (6d) one finds that 
they are proportional to the LO Born term: 



M(6c) = Bg, J2 ^ + + 2Cf ^4/6 (4.5) 

and 



fe=i 



(6d) 



(4.6) 



fe=i 



For the other two vertex insertion diagrams we also 



obtain simple expressions: 



(4.4) M(6e) = (B„-^£n3Br^+2i3r'^ + Q^«(l+/3') 



and 



(6f) 



fe=l 



-16k} + 4iT''T''mv{p2^3u{Pi)u{P3MP4} x 



fe=i 



(4.7) 



= 3iBgg e'^jsf ^ (8mVs - 1) + 2Cf ^m' 



fe=i 



-fcl6(5mVs - 1)} + UT^T^mvipi^suipi) x 



{P3)v{Pi) J2 ^'{5f' (8mVs + 1) - 25, 



(fe-i) 

5 



fc=l 



15 



/(s2/32))/(2/32 



M(12toV'S- 1)} 



number of universal relations among the various coeffi- 
(4.8) cient functions h" valid for any value of j: 



Turning to the two box diagrams (6a) and (6b) we 
note that extensive Dirac algebra manipulations lead to 
rather compact expressions for the amplitudes. We have 
expanded the box diagrams in terms of seven indepen- 
dent Dirac structures, the same set for each of the two 
box graphs. Then every Dirac structure is multiplied by 
the sums of products of a small set of basis functions 
and coefficient functions. Thus, we have the following 
compact expansion for the two box diagrams: 

2 

M = zreoi^e'={5(p2)7^"(Pi)«(P3)7p«(P4)E/i'^'^l"' 
fc=l 

+ «(P2)7"7"7''"(Pi)"(P3)7M7a7.'^(P4) /'"^^^'^ 

+ ■mv{p2'^3u{pi)u{p3)v{p4) ^ ft^h^t^ 

+ mv{p2h''u(j)i)u{p3)-ff,^iv{p4) ^ ft^hf'^ 

+ mv{p2)Yi>3l^'u-{pi)u{p3)lt,luv{pi) ^ ft^hf'^^ . 

(4.9) 

There are seven independent covariants in 1)4. 9|l upon 
using the four mass-shell conditions. We have not at- 
tempted to further reduce the set of seven covariants us- 
ing Fierz-type identities which are anyway valid only in 
n = 4-dimensions. Taking parity and the masslessness 
of the initial quarks into account the number of ampli- 
tudes and thereby the number of independent covariants 
in n = 4 is 2 • 2 • 2 • 2/2 • 2 = 4. However, this counting 
may no longer be true in n ^ 4. 

The sums over i in (|4.9|l run from 1 to 15 in the box 
diagram (6a) . Below we list the color factors and analytic 
functions for the two 4-point functions of 1)4. 9|l . For the 
graph jnt) we get: 



-^col — (-'-im-'- mj )(^ kn-'- nl)f 



(4.10) 



where the first parentheses in 1)4. 10(1 corresponds to the 
summation over color indices of the massless fermion line. 
The basis functions read 



.0) 



.3 -zt4^Vi, h. 



~th 



7 ' 



"lO ^ ''"•8 ' 
"14 



2h. 



-2th 



U) 



(4.12) 



The color factor for the second box graph JBJd) is 

Tcoi = (T"r'')(r''r^). (4.13) 

All basis functions are obtained from those in (|4.11(l by 
the interchange {t <-> u), except for the two additional 
functions (with subscripts 16 and 17), e.g.: 



u), 



Ak) _ p,(fe) 



(4.14) 



Ak) 
J3 



j(k) _ 



.(k) 



Ak) 



m). 



(fe-i) 



J7 — ^4 



(k) 



(k) 









Ak) 

Jll — 


fc. 


Ak) 

J12 


= {k- 




Ak) 

J 13 


= (fc-l)cf"^^(t^u 


Ak) 

J 14 


= (fc- 




Ak) 
J 15 






Ak) 

J 16 




<-> u), 


Ak) p,(fc-i) 

J 17 - ^5 



The last two functions fie and fn appear in the ex- 
pansion (|4.9I) only in two sums where h[^^ and h'f'^ are 
present and, consequently, these sums run from 1 to 17. 

One has further relations for the various coefficient 
functions /i^^ which are similar to those in Eq. 14.12|l . 
In this case they are valid for any given value of j except 
for j = I and j = 4. 



'no — 



where z-u 



U) 



2uh!f''/s, h^f 



-uh 



U) 

7 ' 



"•14 

2m2 - 



2h 



(i) 

7 ' 



'2uh 



7 ' 



(4.15) 



u. In case of j = 1 and j — 4: one has 



4^) = 2«4'Vs, 



4'^ = -w4'\ 



"lO — ""8 ' 



'■12 



'13 - 2u/iilVs, 



'^15 ~ "'^14 • 



(4.16) 



Ak) _ „{k) Ak) _ „{k) 

Ji ~ ^1 ^ h ~ ^5 ^ 



(4.11) 



Ak) _ ^(fe-l) 
/a — ^1 J 



Ak) _ ^(k) Ak) 



C. 



(k~i) 



Ak) _ ^(k) Ak) _ ^(k-l) 



Ak) 
J9 



D 



(k-i) 



Ak) 
Jia 



D 



(fe) 

2 ' 



.(fc) _ ^(fc) 

/ii' = K 



/l^-) = (fc-i)cf-^\ 



/^^)^(fc-l)cr^^ ^{^k-\)D'-^ 

As in the case of the gluon fusion boxes there exist a 



Ak-1) 



Ak) 



Jk-2) 



The coefficient functions ft-p'' are given in Appendix B 
of this paper. However, there exists a partial symmetry 
for these box diagrams, which allows one to express most 
coefficient functions for the box graph |(B|d) through the 
ones of the box graph ^^). In particular, starting from 

the coefficients h" with superscript j > 2, we find the 
following general relations: 



(j)r 



(i)r 



J = 2; (4.17) 
j = 3,5,6. 
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Consequently, for the graph (O)) only the coefficients 

h^P , h']^^ and hff^ are presented in Appendix B. We re- 
iterate that all the one-loop amplitudes of this chapter 
must be multiplied by the common factor lfO|l . 

V. CONCLUSIONS 

We have presented analytic 0{£^) results on the one- 
loop amplitudes for gluon- and light quark-induced 
heavy _c[uark pair production including their absorptive 
parts .30|. These are needed for the calculation of the 
loop-by-loop part of the parton model description of 
NNLO heavy hadron production in hadronic collisions. 
We have not included the finite and divergent pieces in 
our presentation since these were already obtained in an 
earlier publication 20] . The advantage of having the re- 
sults in amplitude form is that one retains the full spin 
information of the partonic subprocess which would be 
of later use when one wants to consider polarization phe- 
nomena in heavy hadron production. As an immediate 
next step we plan to square the one-loop amplitudes and 
to sum over the spins of the external partons. This will 
provide the necessary input for the loop-by-loop part of 
the NNLO parton model description of unpolarized heavy 
hadron or top quark pair production which is presently 
under study at the TEVATRON II and will be studied 
at the upcoming hadron coUider LHC. 
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APPENDIX A 

Here we present the coefficients of the box contri- 
butions for the gluon fusion subprocess appearing in 
Eq. JTTl . 

We define a shorthand notation: 

zi = m? s — , Z2 = s + 2t, 

zt = 2m^ + 1, Zu = 2n? + u, (Al) 

D = m^s — ut. 



First we list coefficients for the abelian type of box 
diagram 

i^l^O, z = 1,2, 3, 4, 13, 
bf} = -lOtzt/D, b'^^l = 2tz.JD, 
6g = szjD, = ~5stf/D, 

6(0)^ = -suipyo, b%\ = s{D + m^s(3'')/D, 
= -2Atzt/D- 

feW = 0, 6« = ~2z^JtD^ 6« = 0, 
b^^l - 2zt/D, b^^l = 2tzt{6D + st(3^)/D\ 
foW ^ -^2zt{2m^D - sf[5^)/D^, 
b^^l = {2m^Z2D + s2tzt/32)/D2, 

= IQm^zt/tD, b'^^l^ = UztiQD + st0^)lD^; 

4!2 = -4m2(2i:)2 + 21^0 + st{D - 2f )l3^)/sfD'^, 

foW ^ Uuzt/sD^, b[^l = im^iD - 2szt)/sD^, 

b'j^]l = Uzt{2m^D - t{D + m'^s)/3^)/D^, 

be.l ^ -8m^{T{sT - 2m^t)D + st^ ztP^)/tD^ , 

b'g l = 4:mh{uz2D - s^zif3^)/sD^, 

bi'l = 2t2/32((2m2 - s)D - mh^P^)/D\ 

= Am'^tl3^{{2s + t)D - s^tff)/D^, 
b[l]^ = 16m'^{{2sT + tzt)D -3m^stzt)/st^D^, 
^14,2 = 4izt(4m2£) - t{D + irri^ s)0^) / D^- 

h^^l = -4(2D2 - 2m2(3m2 - s)D + m*stp^)/ sTD'^ , 
ll^l = -/lzt{2m'^ D"^ - t{2m^D - (3s + 2t)ztD 

- 2mh^Zu))/st^TD^, 

b[]l = 4(s2(to2 - + (m^ + 4s)l3^D 

-tzi(3^ + ztD)/s(3^D^, 
b'-^l = -4i(2(2m^ + 2s~ t)D^ ~ t{Qm^Z2 + s'^0^)D 

+ m^s{QsD + 2tD + stZu)f3^)/sD^, 
b'i}l = ~4:{4:m^TD^ - t{{2m* - f )zt - 2mHz2)D 

+ 2m*sfz2(3'^)/tD^, 
bl^]l = -2{2m^sD^ ~ {2m^v? - s^f(3^ - 2mh^)D 

- 2m'^s'^ziuP'^)/sD^, 

bl^^l = -2t{D^ + (Ws - 2mH + st- 3t^)DP^ 
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+ s2(3i?-m2M)/?4)/i:)3, 

b% = -2zt{tztD + sH^P^ 

- (2w?s'^z^ - Tr?stu + t^)l3^)/D^, 
feW^ = l6m^{s{AD + t(8m2 + ?>s))D0^ 

= Atzt{2{2sT - 3sz2 - 3tu)L> 
+ su(I? + 3m2s)/32)/sL)3; 

43 = -4(s^t-D + + tzt)0^)/s^P^D^, 

= -4(3T/32d + - 2t^ztf}^)/sp^D\ 
ftW = -4tet(6£)2 - 4tTD + sHTp'^)/sD^, 
6g) = 42t((2m^ - i2)£> - 2st^T(f)/D^, 
ftW = -2(2m2(3s + 2t)D'^ + ZsHTD0^ 

+ m^stz2D - 2sH^Zt0^)/sD^, 
ftW = 2t(4m^D2 - s{Qm^s - 6mH + 7st + 2t^)D0^ 

ftW^ = 2i(D2 - (3m2s - lOmH + Ast + 1^)00^ 

+ 2sr'p^)/D'', 
feW^ = 16m^Zt{6tD - s{2sT - f)P'^)/sHp^D'^, 
= -Atzt{2{6m'^s - Sm^t + 7st + 3t'^)D 



bYl = ATu/sD, 

hfl = -2T{2D + st)/stD, 

6g = -4u{D - tzt)/s^(3^D, 

6g = -2((4m2 - 3s)D + stzt)/s^P^D, 

6(2) = 2t'^zti4m'^s + tz2)/sD'^, 

6g = -2(2(m4 - tT)D + t^Zt)/D\ 

6g = {2{mH + szt)D - st'^Zt)/D'^, 

hfl = tiAm^uD + st{Aw?s + tz2)0^)/ sD"^ , 

h^^l^ = t{4m\D + s2(2D - t'^)p^)/sD^, 

&^3^i = lerri^is^TP^ + 2,tuzt) / sH0^ D , 

bfl^ = U^ztiAm^s + 2st + 3t'^)/sD'^; 

6g = -4{D-m\)/sD, 



02I = -2(2(s + T)D + m^st)/stD, 

= -At{m'^Z2 + sZu)/s'^l3'^D, 
6g = 2(£>/32 - tZu)/sP^D, 
hf}^ = -2t(2zj(4,s - t)D + m^sh0^ 

+ 3in^stz2 — st^ Zu) / sD^ , 
hfl = -2{2T^D + t'zu)/D'', 
hf)^ = -(2£)2 - 2mHD + euz2)/D^, 
45 = -ti^D + (6sD - Sm^si + 3sH - 3t^)P^)/D'^, 
b% = -f{4m^D - s\p^)/sD^ 
43,2 = 16(s((m2 + s)D + 2m^st)0^ 

-mh^(2T -u))/shp^D, 
6(^)2 = -4:tzt{10sD - 2rr?sz2 + 3t^u)lsD'^; 

hf^^ = 4{m'^s + TD)/stD, 
bi^l = 2{2{s + T)D + ni^siArr? + t))/stD, 
bf^^ = -4{sD + m^s^ff - mHz2)/s^i:fD, 
bf^^ = -2{2sD + 3w?s^p'^ - 2m^tz2)/s^f3^D, 
bf}^ = 2{4m'^D'^ + 4t(6m'^s + m^t + st)D 

+ st^{5s + 3t)0^ + 2t^zi)/sD^, 
bfl = 2{{2m^T + t^)D - + 2m*te2)/£'^ 

b^^l = ((24m2s - B,mH + 7st)D 

- sH{3m^ - 2t)0^ + 2t^zi)/D'^, 
bfl = (2m^{s - t)D - 2w?sHf3'^ - st^Zu)ID'^, 
bfl = -t{uD - ((8s - 3t)D + m2s(5s + &t))0^)/D'^, 
bfl^ = {2{s'^T + 2mh'^)D - s^{2m^zi - t\)p'^)/sD^, 
brli = 16((16rn* - s'^)D + mHu{^Zt - Az2))/sH0^D, 
b^H^ = Uzt{{lQs - 3t)D + m^s{As + 5t))/sD^; 

bfl = AT{m'^s + D)/stD, 

bfl = 2T{3szt + 2t^)/stD, 

bfl = -4{sDp'^ + mHz2)/s^P^D, 

6g = -2((3sT + + tzt)/s0'D, 

bfl = 2(4m2£>2 - 2t^D + st^(3^ + Am^t'^zi)/sD^ , 

b'^el = 2{2m^TD - 2st^Tp^ - t^Zt)/D^, 

bfl = -t{2tD - st'^p'^ - 4m20l)/£>^ 

ftg = -t{2s'^Tp'^ + 2m^D + stzt)/D'^, 

bfl = i^{Aw?D + s{Am^s - tz2)P'^)/sD^, 

^10,2 = -(2(4m^s + 2mHu + sH)D 



6(3)2 = 4:t'^zt{2m^s - 2st - St^)/sD^; 

b[^l = 0, 41 - 0, 

hfl = -tz2/D, 6(^1 = -{2D + s2/32 + ^^^yjj^ 
b[l, = -24/sP^ b[t, = 0; 

ftg = AT/tD, 6g = 4T/tD, 

= Azt/sDp'', bf} = Azt/sDI3\ 
bfl = A{TD + 2ezt)/D\ bfl = AtTz2/D\ 
b^^} = 2t{D + 2zi)/D\ 
6g = _2i(£) _ 2s^/?2)/£)^ 
6j'oli = -2T(2D + s2/32)/£,2^ 

= -24(m2s/3^ - tZu)/stD0^, 

b% = 12t22*/D2; 

&g = Am^/tD, b^^l = Aw?/tD, 

b'il = 4zu/sD0^ , 6g = 4zJsD(3^ , 

b^il = A{TD + 2mHz2)/D^ 6^ = ^t^zJD^, 

&g = 2{2'm?sz2-tD)/D'^, 

hfl = 2t{D-2sup'^)/D'^, 

bw,2 = '^{{s0^ + '^T)D - m?s'^0^)/D'^, 

^13,2 = -2A{m'^s0^ +tZu)/stD0^, 

b% = -I2tuzt/D'^- 

&M = 0, b'^^l^-2zt/D, 

^3,1 = 0, 4a = ~2{4D - stp'^)/sDp^, 

41 = ^2tzt{iD + tz2)/D^, 

b^Jl = -2ti6TD-st''p^)/D^, 

41 = ./3V(2^t)4a. 
6^o1 = t{8w,^D + 2stu(3^ + tzl)/D^, 
&i3,i(7) = 16(m2s/32 + 2D)/sD0^, 
b^^l, = -4tzt{4D + tZ2)/D^; 

b^^l = -4((2m^s - tz2)D - mH^Z2)/stD'^, 
^2,2 = -4(4£)^ - stD + 2mH'^Z2)/stD'^, 



b^Jl = A{2TD + {2sD + t^u)p^)/sp^D^, 
b^l = -4((3m2z2 - _ 3st)D 

+ 2m^shp^)/s^l3^D^, 
b^l = -4(2D3 - At{2m^ - s)D^ 

+ 2t^u{2s + 3t)D + m^s^t^P^)/sD^, 
41 = M{6m* - 2m^s + 'imH + t^)D 

-2m^st^0^)/D^, 
6^^ = 2(4^tZ)2 _^ 2u{2m'^s - t^)D 

41 = 2t{t{2sT - 12mH - s^)D 

+ s\A{m^ - s)D + i2^„)/32)/s£»3^ 
4o,2 = -2(12m^i£)^ - s^(2m^s + 1^)£>/3^ 

- 2mH'^Z2{D + s^f3^))/sD^, 
^13,2 = 16(t(2m^u - S02)£> 

+ s(5m^s£> + t^D - 3mH\)/3'^)/sH0^D^, 
6^41 = 4:tzt{2s{2zt - 3s)D - 2,t^{s^0^ - tz2))/sD^ 

b^^l = -4{2ztuD - mHz2)/sTD^, 
b^Jl = ~A{ATD^ - t{Am^u - st)D 

+ 2mH^Z2)/stTD^, 
6g = 4((2T + s)D + (m^s^ + fu)P^)/sP^D^, 
41 = -4((m2s + 6m^t + s^)D 

+ 2sHUI3^)/s^P^D\ 
b^^l = At{2{Am^ - u)D^ 

- t{s^ - 2t^)D + m^sHul3^)/sD^, 

b^Jl = 4i((6m^ - 2mH - t^)D + 2m^ stu/S'^) / , 
b^Jl = -2t{AD'^ - tztD 

- {{Am^s - 3s^ + 3i^)£) - m^suz2)l3^)/D^, 
6('^)_3 = 2{8m'^uD^ + m^stz2D 

+ sH{m^ - t)D0^ - 2m^sHuz20^)lsD^, 
5(^)3 = 16(2tM(3m^ - s)D 

+ rr?{AsD + Ztn^)sfp-)l sHfi'^D'^ , 
6(^3 = Uzt(i%m^s - 2st + 6f2)D 

-^sH^"^ +3st^Zt)lsD'^; 

41 = -AT{2sD - t^Z2)lstD^, 

41 = -8T(s£l + t^Z2)lsW'^, 

6g = -4(2(3rn2 + 2t)D + e0')ls0^D'^, 

41 = -4(m2(s - 2€)D - 2st' f3^) / s"" p"" , 
41 = -4(2£)3 _ 8^2^£,2 



- 2t^{2s + 3t)D + sh^T 13^) / sD'^ , 

h^^l = U({Qm^T - 2sT + t^)D - 2st^T(i^) /D^, 
41 = 2{l2m^D'^ - {Am^sT - 3t^Z2)D + st^/3^)/D^, 
b^Jl = 2t{iD'^ - 2m^ztD 

+ (3(m2s + f)D + st^zt)l3'^)/D^, 
=: 2(2(2m\ + szs)!?^ - 2m^t^Z2D 

+ s^tztiD + 2f)f3'^)/sD^, 
b'^^^\^ = 16i2mH{2s + M)D 

+ 'is{D^ + m'^t^)P^)/s^tp^D^, 
b[2^ = 4tzt{2{Am^s + 2st + 3t^)D + 3sfzt)/sD^; 

b'^ll = -iTuz2/sD^, 

b'^Jl = -4(2D2 + t{D + 2mh)z2)/stD^, 

b'^^l = -4(2m^D - t^up^)/sf3^D^, 

b^^l = A{{sT - Qmh)D - 2m'^ s^t(3^) / , 

b^^l = UztiiD^ + m^stz2)/sD^, 

blJl = -4t(2D2 _ (Sto^T + tzt)D + 2m'^st^P^)/D^, 

bill = ~2i(4(2m2 - 5)1)2 ^ 2mhz2D 

- m^s^{2D + tz2)P^)/sD^, 

fo^^^g = -2(2(2to2 - s)z2£'^ + 2m^stztD 

- s^f{D + 2m^Z2)0^)/sD^, 

b^^l^ = lQ{2mh'^D + s(2D2 - 3mh^u)P^)/shp^D^, 

&^4^5 = AtztiSD'^ - tz2£) + 3m^stz2)/sD^. 

Next we list the coefficients for tlie nonabelian box 
diagram |j2t2): 

bf\ = 0, bf\ = 2z,/t\ bf} = 0, 

C = 0, bl°} = -s{4T + s)/D, 

bfl = -16tVi?, 6^"^ = -2(szt + 4t2)/D, 

fefo^i = s{s - 2t)/D, b^°l^ = siszt + At^)/D, 

bill = -^Qrnyt^; 

b['} = 0, bi'} = 2zt/D, 
&W=0, bil^2{m^s + f)/sD, 
= {2{sT - AmH)D - t^zl)/D^, 
h^^l = 2t^{8D ~ szt)/D^, 
bi]l = 2t{ATD + t^Z2)/D^, 
bill = (D^ - ^t'D + sfz2)/D^, 
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6^ = -ATuzt/tD'^, 

^2.1 = -4m2(£) + 2tzt)/tD^, b^^l = -At^u/sD^, 

6^^^ = 4m2(-i:) + 2stl3^)/s(3^D^, 

6^ = 2((4m45 - 2m2t2 + sf)D + m^sh^P^)/D\ 

be.l = Am^{-2m^D^ - shD0^ + s^f (3'^) / sl3^ , 

b^il = ^t^i^m^sT - fu)/D^, bl^l = 8m^fzi/D^, 

bi3,2 = ^Gm^iD + 3m^s(3^)/sP^D^, 

b[W = Af{TD + 3m^zi)/D^; 

b'^^l = Am^Zt{2D + mH)/tTD'^, 

b^^l = -4((2r(5m2s + AmH - f ) + 3t^Z2)D 

- 2m'^t^z^)/stTD^, 
b^^} = AtuysD', 

b[]l = 4((u(8to2 - 3s) + sT)D - 2m'^ uP"^) / 0^ D"^ , 
b^^l = 2(4tL>2 + ((lOm^s + f){2m^ - s) - 2sh)D 

b^ll = 2(2TO2(3zt + S + t)D^ - s^{2m^s - f)D0^ 

-2mh^l3^)lsl3''D\ 
b^^l = At^{{3s{2w? - s) + 2<Z2)i:' 

- rn^s^p'^ + stz2) / sD'^ , 
b^^l = AtztiW ~ 2w?su)lD^, 

b^g^g = -16m2(4(i(5m2 + 2t) - s^P^)D 

+ s{2m'^su - s^zt + t^)(3'^)/s'^tp^D^, 

foW^ = Af{8D'^ - s(7s + 5t)D 

+ 3m^s{~s^l3^ +tz2))/sD^; 

b^^l = ATzt/D'', b^^l = 8tTz2/sD^, 
bi'} = AtysD^ 

6^^] = -4((3m2 + 2<)i:i + 2t^l3^)/sP^D^, 
5^ = -2(4(3m2 + 2t)D^ 

+ t^{2m^ -s- At)D + st*(3'^)/D^, 
bl^l = 2(2(2to* - szt)D^ - s^{2ni^T + tzt)D0^ 

- 2st^{m'^zi ~ t^zt)P^)/sP^D^, 
6^ = -Af{6D^ + stD - t^Z2)/sD^, 
4!1 - -4i2((2r + t)D + 2fzt)/D'', 

6^3^_4 = 16m2(2(2T + s + 2<)i:i - 3sfp^)/s^f3^D^, 
b[W = -4t2(12L)2 _ 2tuD - 3t^Z2)/sD^; 
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"i5 — 



"i2 ' 



(A3) 



0, 



^(2) 
2,1 

f.(2) 



-2T(3s + 4i)/si:i, 



(2) 



^(2) 
8,1 

(2) 

10,: 

(2) 

12, : 

(2) 

13, : 



f,(2) 
^4,2 

^(2) 



^(2) 



^(2) 

^12 

^(2 
''13 



(3) 



b 

b^' 



(3) 



b'i 
b^' 
b 

b 



12,1 

(3) 
13,1 



J,(3) 
''1,2 



0, 6^^i = 2t2)/32_ 4^2 ^^)/^^2_j)^ 

(-2(2771^ - s)D^ + 2s{m^s + tT)D0^ 

2t^{8m^s + 5st + At^)lsD'^, 

-2t{2sTzt + t^)/D^, 
= -{2m^s^T - AfD + t^Z2)/D'^, 
= -t^{2sD + t{Sm^s + bst + At^))/D'^, 
= -lQ{sTu(3^ + 2m^Z2)/s^P^D; 

0, 

2(2(m2s - 2t'^)D - mH'^{3s + At))/st^D, 
0, 

{tz2 - (Sm^s - 2st - t^)l3'^)/s0^D, 
{2mHz2D - s{szt - AtT)D(i^ 

- mhHz20^ - 2m^shp^)/sp^D\ 
-2i2(8(s - t)D - tu{3s + At))/sD'^, 
-2{D^ + ht^D + m^st{Am^ + 

= -{2m'^s^T - At^D - t^u{Zs + 2t))/D'^, 

= -t^{2{s + At)D + tu{3s + At))/D'^, 

= ^{-2sD^0^ + st{2mH{3s + 2t) - 2,uD)0^ 

- 2mH^Z2)/sH^0^D; 

0, 

-2(2(m2s - 2tzt)D + mH'^{5s + At))/st'^D, 
0, 

-(Sm^L) + s{7m^s + 3t'^)0^ + Zt) / P"^ D , 

+ 3m^z2D - (lOm^s - AmH 
+ 2st - 31^)00^ + m'^st{2s0^ + 3z2)0^) / 0^0"^ , 
2t^{16sD + tu{5s + At))/sD^, 
2(1)2 ^ 5^2^, _^ ^2^2(-5^ At))/D^, 

= -(2£)2 - iQstD - Am^sf - bst^u)/D'^, 
= tiAD"^ - 5stD - m^sti^s + At))/D^, 
= 8(t(16m^ + su)D + ni^s'^{2D - tiAm^ - t))0^ 
+ mH'^{hs + At)z2)/sH'^P'^D; 

0, = 2T{8m^s + 3st + At^)/stD, 



u(3) _ 
"3,2 — 

^(3) _ 
"4,2 — 

^(3) _ 
"6,2 — 



^(3) _ 
"7,2 — 

^(3) _ 
"8,2 — 

- 

"U0,2 — 

"12,2 — 

J3) _ 
"13,2 — 



0, 

{-{D + s{Am^ + 1))0^ + 3szt)/s0^D, 
(12m2£)2 ^ 2t{m'^ - s)z2D 
+ st{Am^zi - stzt)P^)/s/3^D^, 
-2t^{5s + At)/sD^, 
2t^{2D + AmH-t^)/D^, 
-- - {2m'^s^T + 3st^ + 2i^) /D"^ , 
--t{2sztD + t^{5s + At))/D^, 
-16(s((4m2 + t)D + mH^)0^ 
+ QmHz2)/sH0^D; 



"1,1 — 
6^^^ - 

"4,1 — 
"7,1 — 
"l3,l - 



J5) _ 
"1,1 — 

J5) _ 
"4,1 — 

"7,1 — 
(5) _ 



"13,1 — 



0, = 2/t, b^il = 0, 

2/sP\ = (D + sz2P^)IP^D, 

0, bfl = -2st/D, 
-- 8(4m2 + M)/st/32; 

0, b'il = 8T/tD, 6g = 0, 
8zt/s(3^D, b'i} = 2{Z2D + 2stztP'')/(3^D\ 
8t^/D^, bf} = -At{D - 2t^)/D^, 
32{-m^sl3^ + tZu)/st0^D- 



A5) _ 
"1,2 - 

z,(5) _ 
"4,2 - 

J5) . 
"6,2 - 

"7,2 - 

J5) 
"13,2 



"1,1 

"4,1 
"6,1 
"7,1 
"8,1 



J7) _ 

"1,2 — 

"2,2 — 

rSJ) _ 

"3,2 — 

vp) _ 

"4,2 — 



J7) _ 
"5,2 — 



0, 43 = SmVii), 43 = 0, 
8zu/sp^D, 

-222(D-2m2s/32)//32_D2, 
-8t\/D^, = + 2tu)/D'^, 

= -32(7X1^ s 13"^ +tZu)/stf3'^D; 

0, = 2t/D, b^J} = 0, 

~2zt/l3^D, 

-{AP^D^ + sztD + stziP^)/P^D^, 

-2f{AD + st)/D^, 

-2st^/D'^, 6^2i = 16m2z2/s/32£>; 

-ATu/D'^, 

-4(6£)2 ^ ^(^g^ ^ 3y)£, ^ 2t^u)/stD^, 
A{ztD + m\2D + tz2)(3^)/sf3^D\ 
-A{3m'^Z2D + {AmH - 3szt - As^0^)D0^ 
+ 2mhtz2p'')/s''p^D\ 
2{ztD'^ + sTztD0^ + (8m*s - Sm^s^ 
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h^^l = 2{-l2m^Z2D'^ + s^{2w?zi 

- 5m^sz2 - sf )DP^ + m^siUD + s^Z2)D0^ 

fe^^) = At^iAD"^ + 2u(2s + t)D + m^s^t)/sD^, 
bsl = 4t((4TO2s + 2mh + t^)D + 2t^u)/D^, 
b[l\^ = l&{lQmHztD + 2m^s^{2D + m^s)P'^ 

- sHi'im'^s + 7mhzu + 2tu^)0^)/ sh(3'^D'^ , 
^15,2 ^ 4:t^{2i-s^P^ + u{3s - 2t))D + ■ist^u)/sD^; 

h^^l = Am'^{2D + mh)/TD^, 

^2,1 = -4(4Ti:>2 + t{2m^T - 3st - 2t^)D 

+ 2m^t^u)/stTD^, 
b^Jl = iz^iim^D + s{2D - tu)P^)/s^P'^D^, 
b^Jl = 4:{2m^{4:m^Z2 + szt)D - ^.s^im^ ~ s)D0^ 

6^^) = 2(-2i:>^ + 2tztD^ - 2m^{3s + 2t)D^P^ 

+ 2s{4m^ - 3s)D^/3^ - st{s^ + AmH)D(3'^ 
+ 2m^st{&m^ - t)D0^ - rr?s^tzl0^ 

- m^sh^iSs + 2t)(3'^)/sP^D^, 

b^Jl = 2{3m^Z2D^ + m^{8m^ + 3s + m)D'^0^ 

+ sh{lQm^ + t)D(j'^ + 2m^s^{zt + 2u)D(3'' 

- 2m^sHzl0^ - 2m^s^t^{3s + 2t)(3'^)/ sfi^D^ , 

h^^l = At^{2{2m^s -s^ + t^)D + shU)/sD^, 
b'^Jl = 4t{{4m'^s + tzt)D - 2m^stu)/D^, 

bi3.3 = I6i2t{8m'^zt + s^Zu)D + lQm^shD0^ 

+ m^s^{AD + 3st)l3^ - 3m^s'^t^ztP^)/shp'^D^, 
b[23 = 4f{2{-s^/3^ + u{s - 2t))D - 3stu^)/sD^; 

^2.1 = -ST{{2s + t)D + st^)/stD^, 
bi^l = -Aztiim^D - st^f3^)/s^f3^D\ 
b[^l = bi^l - 4{3sD + 2tzi)P^)/sP^D^, 
= b'-Jl - 2(4z„i^2 _ 2s2(,„2 ^ 2t)DP^ 

CI = CI - 4(2(2m2z2 + sT)D^ 

- shim's - 3t^)Df3' + sH^(3^)ls(3^D\ 
= 4^^(2(2771^5 + t^)D + shT)/sD^, 

b^Jl = 4t((2sr + tzt)D - 2t^)/D^, 



b^^l^ = bXl^ - 16((47772 - 3s)D + 3tz2Z^)/s0^D^, 
h^^l^ = Af{2{Arr?s + tz2)D - 3st^)/sD^; 

b^;"} = -ATu/D\ 

b^Jl ^ -4:{2D^ + t{2T + s)D + 2t^u)/stD^, 
b'^Jl = 4{ztD - m^iszt - 4D)/3^)/sP^D^, 
b'-ll = 4(3(2777^* - sT)D - 2m^Z2DP^ 

~2st^Zuf3'')/s^l3^D^, 
b'-JJ^ = 2(2(16777^ - 8777^^^ + sh)D^ 

+ 2m^s{3sD - 2m^tz2)D0^ 

+ 277725^(37772 + At)DI3^ - m^st^zl0^)/s^(3^D^, 

b^} = 2{2{hjr?zt - sz2)D^ + 2m^{6s + t)D^P^ 
+ Am^t^Z2DI3^ + st{-im^u - st)DP^ 
+ 2m^s{D^ - st^Z2)(3'^)lsl3^D^, 

b^^l = 4t^{4D^ ~ tz2D + sfu)/sD^, 

b^Jl = U{2D^ + tztD + 2t^u)/D^, 

^Ss = -16(tz„£' + mh{tzt + sz2)f3^ 

+ u{2TD + t^u)P'^)/stp'^D^, 
6^^,. = 4:t^{8D^ - 2tz2D + 3st^u)/sD^. 

Finally, the coefficients for the crossed box |(2t4) are: 

b^l - 0, 4°1 = --*/^^ CI = 0, 

b'^ll = 7t{2D + tu)/sD, 

b^e l ^t{2D + s^ -tu)/sD, 

^10,1 = (2(^ - u)D - t{s^ - tu))/sD, 

bfl-^ = {2{t^ + u^)D + tu{s^ - tu))/sD, 

b[]l = 0, 4!i = 2m^zt/tD, 

foW ^ 0, b[]l = 2Z„(L> - 777 27i) /u^D, 

feW = t{2D + tu){tzu - 6D)/sD^, 

CI = (2(2777^ - t)D^ + 2st^D - t^u^zt)/sD'^, 

= -(2z„d2 - 2stuD + tu^zt)/sD'^, 
b[l\ = -t{4:D^ - 2t^D + t^uZu)/sD^, 

= t{4:D^ + 2tuD - tu^Zu)/sD'^, 
foW^ = lQm'^{-tD + 77727/(t - u))/tu'^D; 

b[]l = 4T7iz„/s£|2, 

^ 4^2(21)2 + t{s - 2T)D - 2t'^uzt)/sfD^, 
bi^]l = AfUz^/suD^, 
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bXl = 4m^((^t - u)D + 2m^u{t - u))/suD^, 
= 2t{2m'^{2mh'^ - 2m^u^ - t\)D 

+ t^u^Zt)/.s^D^ 
feW = -4m2(2(2t + u)D^ + t^u{2t - u)D 

6^^) = 4m'^{2D^ - tu^{2t - u)D - t^u^{t - u))/s^D^, 
&W = -2fu^(2m^D - fzu)/s^D^, 
6^o^2 = 'imH^u{tD + m^s{t ~ u))/s^D^, 
= -I6m'^{2{mh'^ - uT{2t + u))D 

6^2 = -2t{2D + tu){4m^s{t - u) + tu^ Zt) / ; 

ftW = -4{mH'^D + u{m\ - t^U)zt)/stTD'^, 
6^1 = 4(2T(i + T)£)2 + t(2m^s - m^Zt + t^u)D 

+ 2mHTu^Zt)/st^TD^, 
43 = -^tUzu/sD-", 
43 ^ AU{zuD - 2u'^zt)/suD'^, 
ftW = 2t{-l2D^ + 4M(2m^ - i)£>^ 

- u^{f)uzt + st)D - tu^zt)/s^D^, 
ftW = 2{-2{4:mh + 2m,^u + 3t^)D^ - 2t^uD'^ 

+ fu'^(tzt - Am^u)D + 2m^t^u^{t - u))/sHD^, 
b^^l = 2(2(2m2 - 3u)D^ - 2tu^D^ 

+ v^{tzt - Am^u)D + 2mhu^{t - u))/s^D^, 
43 = '2.t^u{<oD'^ - u(t + 2U)D + u^zt)ls^D^, 
6io3 = '^tu{2tD'^ - u{tzt - Am'^u)D 

-2mhu^{t- u))/s'^D^, 
6^3 = -16m^{AuD^ + 2m'^{f + u'^)D 

+ 3tu^zt)/sfuD'^, 
= -2i(2£) + tu){12D^ - 2u{zt + 2U)D 

+ iu^Zt)/s'^D^- 

6^2) = -ATu/sD, b^^l = 2T{2D - tu)/stD, 

6g) = 4{mH - D)/sD, 

bf} = -2U{2D + t^)/suD, 

bfl = -t'{2D^ + IQw^sD + 3t2u2)/s2D2, 

6g = + D)(2m*s2 - t^u^) - 2t^D)/s'^D'^, 

6(2) = (-2(2771^5 + tu)D'^ 

+ 2su{2m\ - st)D - t\^)/s'^D'^, 
6(2) = t\{8D + 2,tu)/s^D'^, 



6fo!i = i(4s£>' + 2i(s2 - 2u^)D + t\^)/s^D\ 
6gi = i-^mh^t -u)- 2t\^)D + t\^)/s^D\ 
6gi = t\2{2mh^ + 1\ + u^)D - t\^)/s^D\ 
6^5^^^ = -8{'imH\ - {2m^{t - u) + st)D)/stuD, 
= -4t'^{2D + w?s){D + m?s)/s'^D'^- 

(A4) 

bfl = A{D-m^u)/sD, 

bfl = 2{2{m^u + 2tT)D + t^U)/st^D, 

6^5 = 4tU/sD, 

bfl = 2{{tZu - 2u^)D - m^u^)/su^D, 

6^2) = -t{2{B>m'^su + IstZu + tu^)D - 3fu^)/s^D'^, 

bfl = {2{2m^s^ + Am^sH + m^stu - st^u - t^)D 

4^2 = (-2(2m^s2 _^ 2m^stu + sTu^ + t^u)D 

bfl = t^u{2{7t + Su)D - 2,tu^)/s^D^, 
bw,2 = (2s(t - u)D^ - 2t^{s - 2u)D - t^u^)/s^D'^, 
b?2.2 = -t{2{2m'^s^ + t^u - tu^)D + tu*)/s^D^, 
45,2 = S{{2m^ su^ ^ 2mH{t^ + u^) - stu^)D 

5(2)2 = _4f(2L) + tu){{7t + 4u)D - tu^)/s^D^; 

6g = 0, 

6(^) = -2((8m2 + 1 - 2u)D + 5mHu)/stD, 
"3,1 ~ ^' 

5(2) = -2{{tZu - 2u{3m^ + u))D + 5m\^)/su'^D, 
b^^} = t{2D + tu){2{7t + 8u)D + btu^) / D'^ , 
bfl = t{2{stzt + 3m^su + 2t{t^ + u^))D 

b^^^-^ = t^{2{'im^s + 5su - t^)D - 5tu^)/s^D^, 
6^)^ = -(4(f2 + u2)£)2 _ 2tu{st - 5u'^)D 

6^5^i = -16{{s{mH - u2) + 4m2u(i - u))D 
- ^m^tu^)/stv?D; 

bfl = 0, b'il = -2T{SD + 5tu)/stD, 

bf}^ = -2{{2mH - 6m^u + Stu)D - bmhu^) / su^ D , 
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,,(5) 

"10,: 
°16 ,: 



(6) 



-t^ {2D + tu) {2D + 5tu) /s^D'^, 
-t^{2{m'^s -s^ + Atu)D + U^v?)Is^D'^, 
= t2(2(4m^s - + tu)D + U^v?)Is^D'^, 
= -(2(2m2s(f2 + _ g2^y _ ^2^2)£, 

= -16{{4m'^u{t -u) + stU)D 
+ 5m'^t'^u^)/stu^D; 



"9,2 — 



= 0, 



i = 1,3,5,9,16, 



-lA, 



.(4) 
'4.1 



t{t-u)/D, 



1/u, 

.(4) 
"15,1 



&W = t(u _ t)/D, 
-4{t - u)/tu; 



b 



-AT/tD, 6g = -AT/tD, 

Am^/uD, b^^l = Am^/uD, 

-4(£)2 ^ rn^st^)/sD'^, = 2^Tu/£l^ 

4(£>2 - m^stu)/sD'^, = 4^\/s£)^ 

= -2tTu/D'^, hfl-^ = -2A{2mH + D)/tuD, 
= -6t'^{2D + tu)/sD'^; 

0, bi'} = -SmVti?, b^^} = 0, 
gJ//?!/), 4*?^ = Atu{2D + tu)/sD^, 



^15, 



2i(2m2?i + i:>)/i:>^ 
= 32(2m2u + D)/tuD- 



.(6) 
'lO.l 



"6,1' 



(7) 



b\ 

bzi 
b^ 

b 



5,2 



"6,2 — 



0, -2,1 

2{m^u- D)/uD, 

t{4D -tu){2D + tu)/sD'^, 

t{6D^ - 2t^D -t^u^)/sD^, 
^t^{2D + u^)/sD^, 
= ?>{D -2rn?u)/uD] 

-ATv?/sD'^, 

4((5m^s - Tu)D + 2tTu^) / stD'^ , 
-At'^U/sD^, 

-A{{2m^s -Tu- Au^)D + 2mHu^)/suD^ 
2t{-8D^ + Au{t + 2u)D'^ 
+ tu^{7t + Au)D + t^u^)/s^D^, 
2t{-QD^ + u^{Am'^s + 5t'^)D 
+ 2t^u^)/s^D^ 



d7) 

-'10,2 
''15,2 



"16,2 

z,(7) _ 

"1,3 - 



,(7) 
2,3 

.(7) 
3,3 



"5,3 



"6,3 



^(7) _ 
"9,3 — 

^(7) _ 
"10,3 - 

"15,3 

"16,3 

"1,4 - 



2t^u{AD'^ - u{3t + Au)D - t\^)/s'^D^, 
■- -2tu'^{{Am'^s + 5t^)D + 2t\)/s'^D^, 
■- 16{2{t - 2u)D'^ + 2u'^{w? - 3t)D 
+ 3m'^t^u^)/stuD'^, 

■- -2t{2D + to)(8£)2 + 8suD - Zt^u^)/ s^D^; 

-A{D^ - m^uD + m^Tu^)/sTD^, 
4:{3TD^ + m^{st - Tu)D + 2m^tTu^)/stTD^, 
AtuU/sD'^, b'-^l = AU{D - 2^M)/sL»^ 
2t{-8D^ - 2u{t - u)D^ 
+ tu^{t - u)D - stu^) / D'^ , 
2t{-6D^ - 2u{t - u)D'^ 
+ tu^{t + 2u)D - 2m^stu^)/s^D^, 
2t\{AD'^ - u{t + 2u)D + tu^)/s^D^, 
= 2tu{{2m^s{t -u)- 3t^u)D + 27X1^ stu^) / D^ , 
= l&w?{2{t + 2u)D - 3tu^)/stD'^, 
= -2t{2D + tu){8D'^ + 2suD + 3tu^) / D^ ; 



(7) 
'2,4 

(7) 
'3,4 

(7) 
'4,4 

(7) 
'5,4 



"6,4 



6^^^ - 

"9,4 ^ 
- 

"10,4 - 
"15,4 
"16,4 



(7) 

1,5 

(7) 
'2,5 

(7) 



"3,5 

u(7) _ 
"4,5 — 

6^'' - 

"5,5 — 

"6,5 — 



J7) _ 
"9,5 — 



AtTu/sD'^, 
4T((3s - u)D - 2t^u)l8tD'^, 
-A{D^ - mhD + mH'^U)/sUD^, 
4((m^ - 3tZu)D - 2mHu)/sUD^, 
-2t((8m'*s2 + 6m^st{t - u) + t^u^)D 

-2t{6D^ + 2u{t - u)D'^ 
+ t^u{At + 3u)D + 2t'^u'^)/s^D^, 
2t^u{AD^ + t{2t + 3u)D + t^u)/s^D^, 
= 2tu{2{t - u)D^ + t^{At + 3u)D + 2t^u)/s^D^, 
= -16(41)2 + 2t(m2 - 2t)D + 3mH^) / stD"^ , 
= -2t{2D + tu){8D^ - 6stD + 3t^u)/s^D^- 

= -4Tw2/sZ)2, 

= 4(£)2 _ rn'^uD + 2tTu^)/stD^, 
= -At^U/sD^, 

= 4(D2 - mhD + 2t^uU)/suD'^, 
-2t{8D^ + AtuD^ - 2t^u^D - m'^sfu^)/s^D^, 
2t{-6D^ + 2{f + u^)D^ + fu^D 
+ 2t\^)/s^D\ 

2t\{AD'^ + tuD - t^u'^)/s^D^, 
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5,(7) 
^10,5 



APPENDIX B 



This Appendix contains the coefficients for the one- 
loop corrections to the subprocess qq QQ. As regards 
the box diagram Fig. Et. we obtain the following coeffi- 
cients h defined in Eq. 14.911 : 



-2T{2/st-l/D), 



hf^ =2(1 + tzt/f3^D)/s, 

- {tzt-sTzi/D + tzt/p^)/D, 



JO) 



-2tT{l + st/D)/D, 



h);" = -t{s^T/D-2t)/D, 
4"^ = {m^s + 2e + st^/D)lD., 
h^^} = 16m\T/t - 2tzt/s^P^)/D; 



h[^'' = -8T/sD, = 8{t + m^Z2/s/3^)/sD, 



h^P = -4ztityD - (1 + l/f3^)/2)/D, 



8t^T/D^ 



,(1) 



.(1) _ /i„+'r/n2 ^(1) 



At{2-t^/D)/D, 



h);' = istT/D^, hYi = -QAm^zt/s^fi^D- 



h[^^=Zt/tD, h^^'' = -l/D, 



/if^ = s(l - st(3^/D)/2D, 
hf^ = -tzi/D'^, 4^' = stz2/2D^, 
^''^^ - szi/2D^, /ifi^ = -8mytD; 



(Bl) 



hf^ = 0, = 0> 

,(3) 



4'^ = szt/4D, hl,'> = ty2D, 
hf^ = st/2D, 4^^ = st/AD, /iJi^ = 0; 



h[^'> = 4T/tD, hi^'> = Azt/sP^D, 
4"*' 2{stzt/D + 2m^Z2/s0^)/D, 
hf^=hf\ h't^=2se/D\ 
hf^ = hf\ hf^ = -16(mVt - z^/s(3^)/D- 



-2/D, 4 = -2z2/sf3^D, 
hf = s{z^lD - Z2lsP^)/D, 

hf^ = 4^' = -l^Zu/s(3'^D] 



.(6) _ f,(5) 



The values for the other coefficient functions with 
« = 3, 5, 9, 10, 12 — 14 and arbitrary j are not written out. 
They can be inferred from the relations presented in the 
Eq. 

The nontrivial coefficients for the second box diagram 
^p) are: 



2{T/D + 2{s + U)/su), 



4°^ = ~2{s/D + 1/s + (2 - uZu/D)/s(3^), 



4°^ = 1 - (Sm^s + 2m2M - s^)/D + ut^{t - u)/D'^ 

?2 



-(2-uz„/i?)//3^ 
4°^ = 2u{m^st/D + to2 - 2u)/D, 
-u{As + 2u- st'^/D)/D, 



(B2) 



,(0) 



4°^ = -2 + s(m2 - 2u)/i:> + m^sH/D^, 
h^°^ = -16to2(.s + [/ + 2tuz^/s^P^)/uD; 

The values for the other coefficient functions h^^ with 
1 = 3,5, 9, 10, 12 — 15 are not spelled out. Again they can 
be inferred from the relations Eq. I|4.15|l . 
Next we write 



,(1) 



4{2mH/u - Z2u)/sD, 



h\'> = 2z2„(l + llP')/sD, 
4^^ -2{2m^s^l3^ + 2utzu + sD)/D^, 
4'^ = 2(zi„(2to2 _ + 2m^Z2u/sP^)/D, 
4^^ = 4u(to2s + uZu)/D'^, 
h^^^ = 2{-2ut^/D + S + 4u) /£», 

4^^ = 2s(m2s-|-uz„)/i:)2^ 
h[\^ = 16m2(3/M - 4z„/s2/32)/£,^ 
h^Y^ = 2u{2uz2u - s'(l + 2/32))/i^^ 



(B3) 



,(4) 



-4/7/wD, /i^^^ = -4z„/s/32l?, 
4"' = -2z2«(m2./i? - IIP^)ID, 
hf^ = -2{suZu/D + 2m^Z2u/sp^)/D, 



,(4) 



= -AuyD^, hf =2sutlD'^, 
4^^ = -2su'^lD'^, h^^^ = lemVuD, 
/lit'' = -2suz2u/D'^-, 
ft^t^ = 4/i?, ht-) = Az2u/sl3-'D. 

The remaining coefficient functions h'f\j = 1,4 with 
i = 5, 9, 10, 12, 13, 15 can be obtained from the relations 
Eq. l|rTH|l . 
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